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ABSTRACT 

The nearby blue compact dwarf, NGC 2915, has its stellar disc embedded in a large, 
extended (~ 22 £?-band scale-lengths) H I disc. New high-resolution H I synthesis obser- 
vations of NGC 2915 have been obtained with the Australia Telescope Compact Array. 
These observations provide evidence of extremely complex H I kinematics within the 
immediate vicinity of the galaxy's star-forming core. We identify and quantify double- 
peaked H I line profiles near the centre of the galaxy and show that the H I energetics 
can be accounted for by the mechanical energy output of the central high-mass stellar 
population within time-scales of 10 6 — 10 T yr. Full three-dimensional models of the H I 
data cube are generated and compared to the observations to test various physical sce- 
narios associated with the high-mass star- forming core of NGC 2915. Purely circular 
H I kinematics are ruled out together with the possibility of a high- velocity-dispersion 
inter-stellar medium at inner radii. Radial velocities of ~ 30 km s _1 are required to 
describe the central- most H I kinematics of the system. Our results lend themselves to 
the simple physical scenario in which the young stellar core of the galaxy expels the 
gas outwards from the centre of the disc, thereby creating a central Hi under-density. 
These kinematics should be thought of as being linked to a central H I outflow rather 
than a large-scale galactic blow-out or wind. 

Key words: galaxies - dwarf, ISM, kinematics and dynamics, evolution 



1 INTRODUCTION 

Quantitatively understanding galaxy evolution is one of the 
major challenges facing modern astrophysics. Star forma- 
tion (SF) activity, along with merger activity and other en- 
vironmental dependencies, largely dictates the evolutionary 
path of a galaxy. Stellar feedback plays a particularly sig- 
nificant role in determining the distribution and kinematics 
of a galaxy's inter-stellar medium (ISM) on local and global 
length-scales. In the s o -called "galacti c fountain scenario" 
([Shapiro fe Field! 1 19761 : iBregmanl [l98Qh , the ISM is heated 
and pushed out of the galaxy by stellar winds and super- 
nova explosions. The ejected material moves into the halo 
of the galaxy, expands, cools and falls back onto the disc. 
Within the gaseous disc, stellar winds can dictate the ISM 
kinematics by setting up expanding gas shells and holes, 
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typically with diameters of 5 - 0.5 kpc (e.g. IWalter &; Brinks! 
I1999I : iBrinks fe Baiaialll986l : lKamphuidll99i for the cases 
of IC 2574; M31; M101 and NGC 6946 respectively). In the 
case that the ejected gas is accelerated to velocities larger 
than the escape velocity of the galaxy, significant fractions 
of the ISM will escape into the inter-galactic medium (IGM) 
(|Mac Low fe Ferrardfl999h . In this sense, stellar winds play 
an important role in the chemical evolution of galaxies and 
the IGM. 

Dwarf galaxies generally have shallow gravitational po- 
tential wells, meaning that the input of kinetic energy into 
the ISM from high-mass star formation and supernova explo- 
sions can significantly affect the gas dynamics. With their 
typically low rotational shear rates and gas surface densi- 
ties, they serve as useful laboratories to study the long-lived 
effects of stellar feedback processes. In this paper we inves- 
tigate the H i dyn amics of the nearby (~ 4.1 ± 0.3 Mpc, 
iMeurer et all l2003h blue compact dwarf NGC 2915 using 
new Hi synthesis observations obtain ed with the A u stralia 
Telescope Compact Array (ATCA). IMeurer et all (p.994) 
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Figure 1. Hi total-intensity map of NGC 2915 (top panel) 
with a zoom-in o f the c entral region (bottom panel), taken 
from Elson et al ] (l2010bh . Contour levels run from 0.08 - 
0.23 mJy beam -1 in steps of 0.02 mJy beam -1 . The intensity 
scale of each panel is specifie d in units of Jy beam -1 . The photo- 
metric centre determined by Elso n et al ] (12010b! ) is marked with 
a cross. The hatched circles in the lower corners represent the half 
power beam width of the synthesised beam. The ellipses have a 
major axis length of 98 arcsec, and represent the i?-band R25 ra- 
dius. The rectangle delimits the plume-like H I feature mentioned 
in Section [3] 



showed the optical appearance of NGC 2915 to be domi- 
nated by two main stellar populations: a blue compact core 
and a more diffuse, older red population. The blue core is 
the locus of ongoing high-mass star formation. This galaxy's 
seemingly typical nature, as suggested by its optical ap- 
pearance, is contrasted by the presence of a large, extended 
(~ 22 ff-b and scale-lengths) H I disc with well-defi ned spiral 
structure (jMeurer et al.lll996l: lElson et alJfeOlObh . The Hi 
total intensity map from lElson et al.l (|2010bf ) is reproduced 
in Fig. [1] with an ellipse that delimits the edge of the stel- 
lar disc. NGC 2915 is also extremely dark-matter-dominated 
(Me urer et al.|[l996h with a total mass to blue light ratio as 
high as Mtot/LB ~ 140 Mq/Lb at the edge of its Hi disc 

(jElson et alJl2010bh. 

Both lMeurer et all (|l996h and lElson etafl (|2010bh pro- 
vide evidence of complex Hi kinematics near the centre of 
the galaxy. The sharp kinks and wiggles seen in the iso- 
velocity contours of their H 1 velocity fields are indicative of 



non-cir cular velocity components within the gas. lElson et al.1 
(2010a) found indirect kinematic evidence for an axisymmet- 
ric radial outflow of 5 — 17 km s _1 near the centre of the 
H 1 disc. Direct evidence of broad and multi-component line 
profiles is provided by position- velocity slices ex tracted from 
the ce ntral regions of the H 1 disc of NGC 2915 ([Elson et al.1 
2010b). These non- Gaussian line profiles are in the imme- 
diate vicinity of the galaxy's star- forming core, and could 
therefore be due to the effects on the ISM of stellar winds 
and supernova explosions associated with the young stellar 
population. In this paper we present the results of a detailed 
study of the kinematics of the central H 1 characteristics of 
the galaxy, with the intention of identifying and quantifying 
the underlying physical processes that lead to the observed 
H 1 properties. A simple physical scenario that links the ef- 
fects of high-mass star formation at the centre of the galaxy 
to the observed H I kinematics is proposed and tested. 

The layout of this paper is as follows. The H i-line data 
set is presented in Section [2] Double-peaked H 1 line profiles 
are identified and quantified in Section [3] with a correspond- 
ing discussion of the H 1 energetics appearing in Section 2] 
Various three-dimensional models of the H i-line data are 
compared to the observations and discussed in Section [5] 
Based on the observed Hi kinematics, the possible fate of 
the central H 1 component is discussed in Section [6] Our 
general conclusions are presented in Section [3 



2 H 1 DATA 



In this work we utilise the H I-line data of lElson et al.1 
(2010b). Using six different antenna configurations of the 
ATCA, nearly 100 hours worth of on-source observations 
were obtained for NGC 2915 between 2006 Octobe r 23 and 
2007 June 2 (project number C1629). lElson et al.1 (|2010bh 
describe the various telescope setups of the observations 
as well as the standard data reduction procedures car- 
ried out to arrive at the calibrated, continuum-subtracted, 
deconvolved Hi data cubes. This analysis makes use of 
their naturally- weighted Hi data cube which has an an- 
gular resolution of 17.0 arcsec x 18.2 arcsec and a chan- 
nel spacing of 3.2 km s _1 . The noise in a line- free chan- 
nel is Gaussian-distributed with a standard deviation of 
a ~ 0.6 mJy beam -1 . 



3 DOUBLE-PEAKED H 1 LINE PROFILES 



lElson et al.1 (|2010bh demonstrated the existence of multi- 
component and broad Hi line profiles within NGC 2915's 
H 1 disc. Many of the profiles within ~ 150 arcsec of the cen- 
tre of the galaxy are observed to be double-peaked. In order 
to classify a particular line profile as being either single- or 
double-component in nature, the reduced \ 2 goodness-of- 
fit statistics of fitted single- and double-component Gaus- 
sians were compared. A line profile with a double Gaussian 
goodness-of-fit statistic, Xdi that was less than 85 per cent of 
the corresponding single Gaussian goodness-of-fit statistic, 
Xs (i.e. Xd < 0.85Xs) was classified as being double-peaked. 
Profiles for which Xd > 0.85Xs were classified as consisting 
of a single component. The Hi data cube was re-gridded to 
a pixel size of 10 arcsec x 10 arcsec in order to increase the 
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level of independence between adjacent line profiles. Three 
niters were used simultaneously when fitting the line pro- 
files with single- and double-component Gaussians: (1) line 
profiles with a fitted peak less than 4cr above the noise were 
rejected; (2) line profiles with a fitted line- width less than 
the channel width were rejected; and (3) fitted profile peaks 
had to be within the velocity range of the data cube. An 
example of a fitted double-peaked profile is shown in Fig. [2] 
A comparison between the reduced \ 2 goodness-of-fit 
statistics of both sets of fits showed that of the 2571 fitted 
profiles, 561 have Xd < 0.85Xs, i.e. they are significantly 
better fit by a double-component Gaussian than by a single- 
component Gaussian. These profiles are shown as black- filled 
squares in Fig. [3] while the profiles for which Xd > 0-85x2 
appear in grey. It is evident that the double-peaked H I pro- 
files are most prevalent within the immediate vicinity of the 
galaxy's star- forming core. Double-peaked profiles also oc- 
cur within the general vicinity of the plume-like H I feature 
seen in the Hi total intensity map (delim ited by a rectan- 
gle in Figs. [Tl and l3|) . lElson et al.1 f|2Q10bl ) show this feature 
within the H I distribution to correspond to asymmetric H I 
line profiles on the receding side of the galaxy. In position- 
velocity slices these Hi profiles appear as an 'Hi beard' that 
is lagging in velocity relative to the main disc. The authors 
estimate the mass of the plume to be ~ 5.6 per cent of 
the total Hi mass. Although kinematically interesting, the 
plume-like Hi feature is not considered any further in this 
work. Also evident from Fig. [3] is the fact that most of the 
double-peaked line profiles near the centre of the galaxy are 
located close to the south-eastern central H I concentration 
seen in the Hi total intensity map. For each of the identified 
double-peaked Hi line profiles within 98 arcsec of the dy- 
namical centre (i.e. within the black ellipse shown in Fig. [3}, 
the absolute difference between the two fitted peak veloci- 
ties was determined. The average velocity separation for the 
profiles is AV — 59 km s _1 , roughly 72 per cent of the 
gala xy's asymptoti c circular velocity of 81.9 =b 1.6 km s _1 
(jElson et al]l2010bh . 



4 CENTRAL ENERGETICS 

Are the effects of stellar winds and/or supernovae explosions 
the cause of the double-peaked H I line profiles observed near 
the centre of NGC 2915? Several authors have linked the 
effects of stellar feedback to radially expanding gas com- 
ponents in nearby galaxies. For example, by carrying out 
a kinematic study of the irregul ar dwarf galaxy NGC 4816 
using Hi and Ha observations, Ivan Evmeren et al.l ([2009) 
detected prominent H I and Ha outflows with expansion ve- 
locities of V exp ~ 25 km s _1 and V exp ~ 30 km s _1 , re- 
spectively. They related these outflows to a fast-expanding 
super giant gas shell, set up by the injecti on of mechanical 
energy from massive stars into the ISM. IWalter &; Brinks! 
(1999), in their study of the holes and shells in the ISM 
of the dwarf galaxy IC 2574, determined the expansion ve- 
locity of the gas around each of 40 identified HII regions. 
They found expansion velocities of V exp ^6 — 25 km s _1 . 
One of their main findings was that the energy requirements 
for "the formation of the holes can be understood in terms 
of the combined effects of stellar winds and multiple super- 
nova explosions of the most massive stars formed during a 
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Figure 2. Example of a double-peaked line profile extracted from 
position CK2000 = 09 h 26 m 12.6 s , 5 2 ooo = -76° 37' 37.8" in the 
re-gridded H I data cube of NGC 2915. This position is marked 
by a black cross in Fig. [3] In the top panel, open circles represent 
the data, the solid grey and black curves represent best-fitting 
single- and double-component Gaussians, respectively. The solid 
grey and black curves in the lower panel represent the residuals for 
the fitted single- and double-component Gaussians, respectively. 

recent phase of active star formation" . lYoung et al.1 (|20Q3h 
used VLA observations of three dwarf galaxies to analyse 
their H I line profile widths and shapes. The authors found a 
trend between the fraction of asymmetric profiles in a galaxy 
and its Ha luminosity, with galaxies having greater star for- 
mation rates also having a greater fraction of asymmetric, 
double-peaked line profiles. They interpret the result as be- 
ing indicative of the star formation activity and its associ- 
ated injection of kinetic energy into the ISM that is respon- 
sible for stirring the surrounding Hi. 

4.1 Gas energetics 

Supposing that the double-peaked H I line profiles observed 
in the central regions of NGC 2915 are due to an expanding 
gas component, could the young stellar core be supplying 
the energy required to drive such an expansion? To begin 
answering this question, an estimate of the ene rgy associated 
with t he expanding gas is required. Following iMeurer et al.l 
(1996), the central energetics are quantified by estimating 
the kinetic energy of an expanding gas component as 

E k = ±M g AV 2 , (1) 

where M g is the mass of the expanding gas and AV is 
the observed velocity separation between the peaks of a 
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Figure 3. Comparison of the quality of single- and double- 
component Gaussian fits to the Hi line profiles. Black-filled 
squares denote the positions at which the reduced x 2 statistic 
of the fitted double Gaussian is less than 85 per cent of the corre- 
sponding statistic for the fitted single Gaussian, i.e. Xd < 0-85x 2 . 
Grey-filled squares represent profiles for which Xd > ^-^Xs- 
White-filled squares represent line profiles that were rejected by 
the fitting filters and which were therefore not parameterised. The 
angular size of each coloured square is 10 arcsec x 10 arcsec, cor- 
responding to « 183 pc x 183 pc. The single contour is at an 
Hi column density level of 1.6 x 10 20 cm -2 , and is arbitrarily se- 
lected to represent the approximate edge of the outer H I disc. The 
ellipse represents the edge of the old stellar disc and has a semi- 
major axis length of 98 arcsec. The rectangle is the same as the 
one shown in Fig.^and delimits the plume-like H I feature seen in 
the H I total intensity map. The cross at a 2 ooo = 09 h 26 m 12.6 s , 
£2000 = —76° 37' 37.8 /x marks the position of the H I profile shown 
in Fig. El 



double-peaked Hi line profile. The observed velocity sep- 
aration is related to the expansion velocity by V exp = §AV. 
For the sake of brevity, the term "splitting" is henceforth 
used to denote the velocity separation between the fit- 
ted peaks of a double-peaked Hi line profile. Considering 
only the split line profiles within 98 arcsec of the centre of 
NGC 2915 (i.e. black-filled resolution elements within the 
ellipse shown in Fig. [3]), the corresponding flux density of 
each profile (obtained from the corresponding resolution el- 
ement in the re-gridded Hi total intensity map) was con- 
verted to a mass estimate, all of which were summed to yield 
M g rsj 2.6 x 10 6 M0 for the mass of expanding gas near the 
galaxy centre. The average splitting of AV = 59 km s _1 for 
the central double-peaked line profiles suggests an expan- 



sion velocity of V e3 



30 km s 



These M g and AV values 



yield Ek ~ 5.6 x 10 44 J for the amount of kinetic energy that 
can be linked to an expa nding central gas component. 

iMeurer et all (|l996h . using M g = 10 8 M and AV = 
40 kms" 1 , estimated E k ~ 4x 10 46 J for NGC 2915 assum- 
ing a distance of 5.1 Mpc. At a distance of 4.1 Mpc, their en- 
ergy estimate reduces to Ek ~ 2.6 x 10 46 J. It is stressed that 
our estimate of M g is a lower limit for the mass of expanding 
gas. Only the masses corresponding to identified split pro- 
files were considered, as opposed to all of the mass within a 
specified central r egion. For this reason , our Ek estimate is 

JI996). 



lower than that of Meurer et al 



rived by lElson et all (|2Q10bh from the global H I profile of 
their high-resolution Hi data c ube. We there f ore tr eat the 
E k ~ 2.6 x 10 46 J estimate of iMeurer et al.1 (|l996h as an 
upper limit for the kinetic energy of the expanding neutral 
ISM. 

4.2 Stellar energetics 

The amount of kinetic energy linked to the expanding gas 
component needs to be compared to the mechanical en- 
ergy output of the central stellar population. Ha emission, 
a tracer of recent and ongoing high-mass star formation, is 
observed near the centre of NGC 2915. Importantly, it only 
occurs at inner radii where a large fraction of the split Hi 
line profi les are observed. An H a image produced using the 
data of Gil de Paz et al. (2003) is shown in the top-right 
panel of Fig. [U The GALEX near ultra-violet image is shown 
in the top-left panel. The location of the Ha emission rela- 
tive to the Hi emission is shown in the bottom panel. The 
Ha emission is clearly contained between the two central H I 
concentrations seen in the Hi total intensity map. The fila- 
mentary structure of the Ha emission is indeed testament to 
the perturbing effects of the high-mass star formation on the 
distribution of the ISM. As an example, the large shell-like 
structure seen in the south-east region of the Ha distribu- 
tion is spatially coincident with the split H I line profiles near 
the centre of the galaxy. 

While the intense ultra-violet radiation energy pro- 
duced by young, massive O- and B-type stars is capable 
of ionising the neutral ISM, it is not expected to play a sig- 
nificant role in redistributing it. Rather, it is the mechanical 
energy injected by stars via stellar winds and supernova ex- 
plosions that is expected to set up shock-heated (T > 10 6 K ) 
gas and expanding structures within the ISM. IChul ([2005) 
neatly demonstrates that the total amount of mechanical 
wind energy output during the entire lifetime of an O-type 
star, including the high-mass-loss-rate Wolf-Rayet phase, is 
~ 4 x 10 43 J. The author further points out that this en- 
ergy is more than a factor of two lower than the mechanical 
explosion energy of a supernova, ~ 10 44 J. In this work we 
therefore ignore the energy deposited into the ISM by stellar 
winds and consider only the mechanical energy produced by 
supernova explosions. In the remainder of this section we 
estimate the time-scales over which the mechanical energy 
deposition from supernovae into the ISM can account for the 
E k ~ 5.6 x 10 44 J and E k ~ 2.6 x 10 46 J estimates presented 
ab ov e for the ISM kinetic energy . 

iLeitherer &; Heckmanl ([19951 ) show that the equilibrium 
ratio of the mechanical luminosity from supernova explo- 
sions to a galaxy's Ha luminosity is L^l^hl 'L hol — 1.8 . This , 
they say, is the case for a Salpeter IMF ([Salpeterl [l955) 
stellar population having formed at a constant rate. Us- 



ing this relation togethe r with t h e Lh o 



1(T 



mass M g 10 8 M used bv lMeurer et all {l996) to estimate 



Ek is almost as large as the total Hi mass of NGC 2915 de- 



measurement of iGil de Paz et al.l (|2QQ3h yields an estimate 
of L s J^l h = 5.7 x 10 32 J s _1 for the total mechani- 
cal luminosity of NGC 2915's supernova population. How- 
ever, not all of this mechanical energy will b e deposited 
i nto the ISM as kinetic energy. For example, iMori et"aD 
(2002) treat only 30 per cent of the mechanical energy 
from supernovae as being converted into ISM kinetic en- 
ergy. Adopting the same conversion factor, we estimate 
L'Jech = 0.3L^ h = 1.7 x 10 32 J s" 1 to be the rate 
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Figure 4. Top left panel: GALEX near ultra-violet image of 
NGC 2915. Top rig ht panel: Hq image of NGC 2915 produced 
using the data from lGil de Paz et al.l (|2003|). Bottom panel: Ha 
image with H I flux density contours overlaid. H I contour levels 
run from 10 - 60 mjy beam -1 in steps of 10 mJy beam -1 and 
from 60 - 230 mJy beam -1 in steps of 20 mJy beam -1 . The Ha 
intensity scale is un-calibrated in the top right and bottom pan- 
els. In all panels the dashed black ellipses, included to facilitate 
length-scale comparisons between the images, have a semi-major 
axis of length 20 arcsec. The solid black ellipse in the bottom 
panel is the same ellipse shown in Figs. [T] and [3] and represents 
the approximate edge of the old stellar disc. The crosses in the 
top panels represen t the photometric centre as determined by 
lElson et al.l (|2010bh . 

at which mechanical energy from supernovae is injected into 
NGC 2915's ISM. At this rate the central Hi energetics es- 
timates of Ek ~ 5.6 x 10 44 to 2.6 x 10 46 J can be accounted 
for within - 1.0 x 10 5 to 4.8 x 10 6 yr. 

A second estimate for the timescale required by 
NGC 2915's supernova population to energise the central HI 
component of the galaxy can be obtained by considering the 
expected total supernova mass and supernova rate. Treating 
a single supernova event as depositing 0.3 x 10 44 J of mechan- 
ical energy into the ISM, roughly 20 - 860 supernovae are 
needed to account for amounts Ek ~ 5.6 x 10 44 - 2.6 x 10 46 J 
of kinetic energy associated with the expanding central H I 
component. As a galaxy produces stars that will later ex- 
plode as supernovae, it builds up stellar mass. A question 
therefore is whether NGC 2915's observed stellar mass is 
consistent with the galaxy having formed 20 - 860 type-II 
supernova progenitors over its stellar assembly lifetime. To 
answer this question we randomly sampled stellar masses 
betw een 0.08 M and 120 M from a Kroupa IMF (jKroupal 
2001). As the masses were sampled, supernova progenitors 
were identified as those stars more massive than 8 M©. 860 
supernovae were found to be accounted for by a total stellar 
mass of only - 6.7 x 10 4 M . These supernovae constituted 
~ 26 per cent of the total stellar mass, with an average 
supernova mass of ~ 19.8 M©. The supernova masses were 
distributed as M^ 12±010 . The remaining 74 per cent of the 



total stellar mass was hence made up by stars less massive 
8 M©. A total stellar mass of 6.7 x 10 4 M© is much less than 
that of the blue stellar population of NGC 2915. Adopt- 
ing a B-band extinction-corrected apparent magnitud e of 
tub = 12.0 (jHelmboldt et al.ll2004l : ISchlegel et al.lll998h to- 
gether with the 4.1 Mpc distance estimate of iMeurer et al.l 
(|2003h . the total B-band luminosity of NGC 2915 is deter- 
mined to be 3.2 x 10 8 L©. Even for a B-band stellar mass-to- 
light ratio as low as 0.1 M©/L©, the total mass of the blue 
stellar population is still at least two orders of magnitude 
larger than that required for ~ 860 type-II supernovae to 
form. 

The supernova rate of a galaxy is often found to be of 
the order of one per cent of the t otal star formation rate 
(jBoom sma 20071) . iKennicuttl (| 19831 ) determined the relation 

SFRtot = 3.42 x 10~ 8 L Ha , (2) 

where SFR to t is the total star formation rate of a galaxy 
in units of M© yr -1 and Lho. its total Ha luminos- 
ity in solar units. Using this relation toge ther with the 
Lhcx — 10 32-5 J s _1 measurement of I Gil de Paz et al.l 
(|2003h yields an estimate of SFR to t = 2.8 x 10~ 2 M© yr -1 
for the total star formation rate of NGC 2915. The su- 
pernova rate is thus expected to be of the order of 
10 -4 M© yr -1 . Considering only type-II supernovae events, 
a lower supernova mass limit for 860 supernovae is Msn = 
860x8 « 6.8xl0 3 M©, since core collapse supernovae events 
only occur for stars at least as massive as 8 M© . The ratio of 
this supernova mass to the supernova rate yields a formation 
timescale of ~ 6.8 x 10 7 yr. The corresponding timescale for 
20 supernovae is ~ 1.6 Myr. 

These estimated time-scales of ~ 10 5 — 10 7 yr match 
the typical main sequence lifetimes of high-mass O- or B- 
type stars. These are exactly the types of stars expected 
to constitute NGC 2915's blue stellar core, thereby impli- 
cating them as a mechanical energy source responsible for 
the anomalous central Hi kinematics. In the sections that 
follow, the central Hi component of NGC 2915 is modelled 
as being radially expanding due to the input of kinetic en- 
ergy from the high-mass stars into the ISM. The kinematic 
models that we produce are not expected to precisely match 
the complex features of the Hi data. Instead, they allow 
us to generally identify the constituent disc properties and 
processes that may be responsible for the observed H I dis- 
tribution and kinematics. 



5 THREE-DIMENSIONAL MODELS 

A radially expanding inner gas component should exhibit 
distinct kinematic signatures as seen in a three-dimensional 
data cube. To check whether such a scenario is consistent 
with the observed Hi kinematics of NGC 2915, various 
model cubes were constructed and compared to the obser- 
vations. The model cubes differ in the radial profiles of some 
of their orientation parameters as well as their kinematics. 
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5.1 Modeling routine 

The model cubes were generated using the task halo- 
MODGAI0, a modified version of the standard gipsy task 
GALMOD. The routine assumes axisymmetry and distributes 
many "Hi clouds" into a set of tilted rings. For each ring 
the inclination, position angle, and rotation velocity must 
be provided; as well as a velocity dispersion, column den- 
sity and scale-height perpendicular to the plane. Unlike the 
standard GALMOD routine, the user can additionally spec- 
ify non-circular velocities parallel and perpendicular to the 
plane for each ring using HALOMODGAL. The routine selects 
random positions at which to distribute Hi clouds within 
each ring. Having placed a cloud, the routine uses the sys- 
temic and rotation velocities of the tilted ring to build a 
velocity profile around the position of the cloud. For this 
purpose the cloud is divided into sub-clouds and each sub- 
cloud assigned a velocity equal to the central velocity of the 
clouc0 plus a random increment extracted from a Gaussian. 
This Gaussian has a mean equal to the central velocity and 
a width equal to the velocity dispersion of the ring. Finally, 
the routine uses the position and velocity of the sub-cloud to 
place it at the corresponding pixel position within the data 
cube. 

5.2 H I disc scale-height 

The HALOMODGAL routine requires an estimate of the H I 
disc scale -height. To obtain an es timate, the method out- 
lined by IWalter fc Brinks! (|l999h was used. They state 
that for a densit y distribution given by p(z, R ) — 
2/o(0, i?) sech(>/zo) |Kellmad[l97ilvan der Kruitlll98lh . the 
scale- height is 

z (R) = —J^ , (3) 

y/2<irGp(0, R) 

where a gas is the velocity dispersion of the gas in quiescent 
regions of the galaxy and p(0, R) is the total volume density 
of the galaxy in the disc plane. Under this approximation, 
the Hi disc scale- height is proportional to the Hi velocity 
dispersion and inversely proportional to the total mass den- 
sity in the disc. Assuming a Gaussian distribution of the H I 
disc in the ^-direction, the la scale-height is approximately 
h = zo/V2. 

A constant value of 10 km s _1 was used for o- gas , chosen 
in accordance with the typical azimuthally-averaged second- 
order moment s of the H i line p rofiles of th e outer H I disc of 
NGC 2915 (jElson et alJl2010bh . Following IWalter fe Brinks] 
(1999) to calculate the average mass density in the plane 
of the gas disc, it was assumed that the total mass of the 
galaxy can be estimated from the last measured point (with 
radius Riast) of the observed rotation curve, Viast- This mass 
was treated as being distributed over a spherical volume 
of radius Riast- The results of the tilted ring m odels fitted 
to the Hi velocity field by Elso n et al . (2010b) were used: 
Riast — 8.7 kpc and Vi as t = 86 km s _1 . These values yield 
an average density of p = 0.005 Mq pc -3 at the radius Ri as t- 
Since this is only the mean density of th e galaxy, and the 
density near the plane is certainly higher, IWalter fe Brinks 

* Kindly provided by Filippo Fraternali. 

§ The sum of the cloud's systemic and rotation velocities. 



(1999) suggest that twice the derived value should be used to 
approximate the mid-plane density. The mid-plane density 
is thus approximated as p(0, R) — 0.01 Mq pc~ 3 . Under 
a minimum disc assumption, El son et al.1 (|2010bh param- 
eterised the dark matter halo of NGC 2915 as a pseudo- 
isothermal sphere. They estimated a central dark matter 
density of po = 0.17 =b 0.03 M pc -3 and a core radius 
of R c = 0.9 ±0.1 kpc. Our estimate for the mid-plane 
dark matter density is therefore roughly a tenth of the core 
density, and matches the halo volume density at a galac- 
tocentric radius of 3.6 kpc. Using these a gas and p(0, R) 
estimates together with Eqn. [3] and h = Z0/V2 yields 
h = 413 pc as the approximate scale-height of the gaseous 
disc NGC 2915. This thick Hi layer is attributed to the 
lower gravitational potential of the system (as compared to 
high-surface-brightness late-type spirals) and the H 1 veloc- 
ity dispersion of a gas ~ 10 km s . 

5.3 Tested models 

Several models were generated to study the Hi kinemat- 
ics of NGC 2915. To create th e se mod els, the tilted ring 
modeling results of lElson et al.l (|2Q10bh were used. Specif- 
ically, the inclination and position angle radial profiles as 
well as the rotation curves of their so-called tilted ring mod- 
els CI and SI were incorporated. These profiles are shown 
in Fig. [5jt, b and c. The two tilted ring models differ mainly 
in their inclination profiles, and very slightly in their rota- 
tion curves. Model CI has an almost constant inclination 
of i ~ 55° for all radii. The inclination profile of model SI 
rises from i ~ 55° to i ~ 75° when moving from the outer 
to the inner disc, beginning at a galactocentric radius of 
~ 240 arcsec. For all model cubes, the systemic velocity was 
always kept constant at 465 km s _1 . The kinematic cen- 
tre was_s^_equaJ__tothe photometric centre as determined 
bvlElson et al.l (l2Q10bh . The H 1 column density profile from 
Elson et al .1 (|2010bh was inclination-corrected using each of 
the inclination radial profiles, thereby yielding the two face- 
on column density radial profiles shown in Fig. [5ji. Each 
of the model cubes was smoothed down to a resolution of 
17.0 arcsec x 18.2 arcsec, the resolution of the actual data 
cube. 

Various position- velocity slices were extracted from the 
model cubes and compared directly to the corresponding 
slices from the Hi data cube. The locations and orienta- 
tions of these slices in the H 1 total intensity map are shown 
in the bottom-most panel of Fig. [6j Each position- velocity 
slice was chosen to investigate a particular component or 
feature of the complex Hi kinematics of NGC 2915. Slice 1 
targets the outer disc and is used to demonstrate its rel- 
atively simple circular kinematics. Slice 2 is placed along 
the kinematic minor axis of the Hi disc in order to probe 
radial and other non-circular gas motions. The general bar- 
like morphology of the inner H 1 disc, including the two H 1 
concentrations near the galaxy's centre, is studied by means 
of slice 3. Finally, slice 4 is selected to further highlight the 
distinct non-circular motions near the centre of the H 1 disc. 

5.3.1 Purely circular models 

Six different model cubes were constructed in total. All of 
the model position- velocity data are shown in Fig.[6]together 
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Figure 5. Radial profiles used to construct the model cubes. Panels (a), (b), (c) and (d) represent inclination, position angle, circular 
velocity and face-on colu mn density radial profiles, respectively. Solid and dashed curves represent the profiles from tilted ring models 
CI and SI of lElson et al.l (|2010br ) respectively. Note that these tilted ring models have the same position angle radial profile, shown as 
a single curve in panel (b). The dotted curve in panel (a) represents the inclination radial profile used to construct model cube Sl+exp 
(Sec.[533}. 



with the observed data. The first two are purely circular 
models, and were made using the radial parameter pro files 
of tilted ring models CI and SI from lElson et all (|2Q10bh . to- 
gether with a constant Hi velocity dispersion of 10 km s _1 . 
These two model cubes, henceforth referred to as cubes CI 
and SI, served as the "base" models. They were manipulated 
in order to construct the rest of the model cubes. Cubes CI 
and SI are shown in columns 1 and 2 of Fig. [6] Slice 1 (row 1) 
demonstrates how these purely circular models are capable 
of reproducing the observed H I kinematics of the outer disc 
where non-circular gas motions are negligible. Slice 1 also 
shows that, by construction, all six models reduce to mod- 
els CI and SI in the outer disc. Slices 2, 3 and 4, however, 
show that simple, purely circular models such as cubes CI 
and SI are unable to reproduce the inner-disc H I kinematics 
of NGC 2915 which are dominated by non-circular velocity 
components. 

5.3.2 Variable dispersion models 

Is a neutral ISM with a high velocity dispersion respon- 
sible for the non- Gaussian line profiles near the centre of 
the galax y? The secon d-order Hi moment map presented 
bv lElson et al.l ([20 10 oh exhibits a sharp increase from outer 
to inner disc, with second-order moments as high as ~ 
25 km s _1 near the centre of the galaxy. Is this sharp in- 
crease in the second-order moments of the H I line profiles 
an accurate representation of the H I velocity dispersions, or 
is it an artefact of the split line profiles near the centre of the 
galaxy? To check whether such an increase in velocity dis- 
persion is consistent with the H I data, the second -order mo- 
ment map radial profile from lElson et al.l (|2010bh was added 
as a velocity dispersion radial profile to model cubes CI and 
SI. The resulting model cubes are henceforth referred to as 
cubes Cl+disp and Sl+disp, and are presented in columns 
3 and 4 of Fig. El Panels 10, 11, 12, 14, 15, 16 of Fig. [6] 
show that although these dispersion models are able to re- 



produce the large spreads in central velocities, they cannot 
match the double-peaked structure of the H I line profiles. A 
central ISM with a high velocity dispersion therefore cannot 
explain the central gas kinematics of NGC 2915. 

5.3.3 Radial flow models 

Since none of the purely circular models were able to match 
the H I kinematics of the inner disc, the final set of models 
included radial velocities. The results of Sec. [3] above show 
that, on average, Hi line profiles are split by ~ 60 km s _1 
within the central ~ 100 arcsec of the galaxy. If this observed 
line splitting is interpreted as being due to an expanding gas 
component, the expansion speed is equal to half of the ob- 
served line splitting. Model cubes Cl+exp and Sl+exp were 
therefore constructed by adding radial expansion velocities 
of Vexp = 30 km s _1 at all radii R ^ 85 arcsec to model 
cubes CI and SI. These models are shown in columns 5 and 
6 of Fig. [6] A distinct trend in the general structure of the 
model cube line profiles is noticeable when moving from col- 
umn 1 to column 6 (left to right) in Fig. [6] columns 1 and 2 
have single-component, relatively low-dispersion Gaussian- 
shaped line profiles. The remaining columns 3-6 have line 
profiles with higher second-order moments, yet columns 3 
and 4 still only contain single-component line profiles. Most 
noticeable in columns 5 and 6 is the double-peaked nature 
of the inner-most line profiles, with the more specific nature 
of the splitting differing between columns 5 and 6. 

The split line profiles seen in the model cubes (pan- 
els 18, 19, 20, 22, 23, 24) are mostly consistent with the 
double-peaked profiles seen in the data (panels 26, 27, 28). 
Clearly, however, it is cube Sl+exp that best matches the 
data. The only significant difference between cubes Cl+exp 
and Sl+exp is the inclination of the inner disc, which was 
found to play a crucial role in determining the double-peaked 
nature of the central line profiles. Various inner disc incli- 
nations in the range 70° - 80°, as suggested by the tilted 
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ring models from Els on et al.l (|2010bl ). were experimented 
with. All such inclinations were found to result in relatively 
poor matches to the observations in terms of split line pro- 
files. Model Sl+exp was finally constructed using a constant 
inner disc inclination of 85° out to a galactocentric radius 
of 140 arcsec. Figure [5ji shows the inclination radial pro- 
file used for cube Sl+exp (dotted curve) together with the 
profiles used for the purely circular models. Without this 
highly-inclined inner disc, the observed splitting of the line 
profiles at inner radii could not be matched. Radial velocities 
of order V exp ~ 30 km s _1 together with a highly-inclined 
inner disc are therefore crucial in terms of the description of 
the central gas kinematics of NGC 2915. 

5.4 Discussion 

The distribution and kinematics of the H I near the centre of 
NGC 2915 are consistent with the simple physical scenario in 
which the young stellar core of the galaxy, with which the Ha 
emission is related, expels the gas outwards from the centre 
of the disc. This process creates a central H I under-density 
and a radially expanding central gas component. Similar 
processes are o bserved in other nearby galaxies. As men- 
tioned in Sec. [4j Ivan Evmeren et all (|2009h detected promi- 
nent Hi and Ha outflows in NGC 4816, with expansion ve- 
locities of Vex P ~ 25 km s _1 and V exp ~ 30 km s _1 , re- 
spectively. Such outflows can significantly affect the overall 
evolution of a galaxy. Star-bursts, with their associated gas 
expulsion processes, are indeed thought to be a possible evo- 
lutionary link between late- and early- type dwarfs. If a burst 
of star formation experienced by a gas-rich late-type dwarf is 
able to remove or at least redistribute a significant mass frac- 
tion of the ISM, star- format ion could be quenched, resulting 
in the galaxy fading to become a red, gas-poor early- type 
dwarf. The redistribution of gas within a galaxy can also 
lead to episodic star formation activity, with the removal and 
subsequent replenishment of H I being linked to correspond- 
ing epochs of diminished and in tensified star formation ac- 
tivity (e.g. IWilliams et aDl2010l ). Indeed, colour- magnitude 
diagrams of resolved stellar populations of nearby galaxies 
provide clear evidence of complex, episodic star f ormation 
histories (e.g. llMonelli et al.ll2003l : IWeisz et alJl2008h . 

Although the high-resolution Hi observations to- 
gether with our modeling results presented above suggest 
NGC 2915 to contain relatively little Hi mass at its centre, 
it could be the case that molecul ar hydrogen (H2) fills the 
central region. iBigiel et al] (|2008h showed that a Schmidt- 
type power law with index N = 1.0 =b 2 relates the observed 
star formation rate surface density to the H2 surface den- 
sity across a samp le of spiral galax ies from The H 1 Nearby 
Galaxy Survey ([Walter et al.ll2008l ). implying that H2 forms 
stars at a constant efficiency. This relation allows the molec- 
ular gas content of a system to be estimated from a measure 
of its star formation rate surface density. Based on the in- 
tegrated star formation rate surface density of the stellar 
core of NGC 2915 to be presented in our upcoming pub- 
lication (Elson et al., 2011), a total molecular gas mass of 
Mh 2 — 5.7^21 x 10 7 M© is estimated to be co-located with 
the stellar disc of NGC 2915. This amount of molecular gas 
is conceivably sufficient to fill the region between the central 
Hi concentrations. 

Unclear is the reason why the inner portion of the Hi 



disc of NGC 2915 should be so highly- inclined relative to the 
outer portion. It may be the case that the disc is severely 
warped. N umerical simu lations of dissipative mergers car- 
ried out bv lBarnes! (2002) have shown that ext remely warped 
gas discs can result from the merging process. Barne s (2002) 
showed further that extremely warped gas disc merger rem- 
nants can result even in the case of unequal mass mergers. 
A warped gas structure can also develop in the case that 
gas from the inter-galactic med ium accretes onto the inne r 
portions of a galaxy's gas disc (|Bournaud fe Combedl2003h . 
The degree of the warping will depend on the ratio of the 
angular momentum of the in-falling gas to that of the disc. 
Although direct evidence for the above-mentioned scenarios 
is lacking in the case of NGC 2915, they cannot be ruled 
out as possible causes of a highly- warped inner Hi disc. A 
good example of a galaxy with an extremely warped disc 
is the "Spindle" galaxy (NGC 2685). Until recently, this 
system was generally regarded as being a two-ringed polar 
galaxy, consisting of two kinematicall y distinct discs ali gned 
almost perpendicular to one another. Ijozsa et~aD (|2009h re- 
cently showed, however, that NGC 2685 forms a coherent, 
extremely warped disc which creates the appearance of two 
mutually inclined regions. 



6 FATE OF THE EXPANDING GAS 

Having established that the central Hi kinematics of 
NGC 2915 are consistent with those of a radially-expanding 
gas component, the eventual fate of the expanding gas can 
be speculated on. Will the central gas eventually stop ex- 
panding and fall back towards the centre, thereby fuelling a 
subsequent burst of star formation, or will it escape into the 
inter-galactic medium? 

Several authors have studied gas o utflo ws in nearby 
dwarf galaxies (see iMeurer et al.l Il992l and iMartinl Il996l 
for NGC 1705 and IZw 18 respectively). In VVI Zw403, 
Pap aderos et a D £[994) detected extended X-ray plumes 
which they interpreted as the outflow of hot gas. 
iMac Low &; Ferraral ([l999) modelled the effects of repeated 
supernova explosions from starbursts in dwarf galaxies on 
the ISM, taking into account the gravitational potential of 
their dark matter halos. They found that the mass ejection 
efficiency is very low in galaxies with a total mass > 10 7 Mq . 
More specifically, they found that in galaxies more massive 
than ~ 10 6 M©, the fraction of mass that is lost in a galac- 
tic blowout is as little as 10 -6 . For refere nce, the dynamica l 
mass for NG C 2915 is ~ 1.5 x 10 1 M© (|Elson et al.ll2010bT ). 
Importantlv JMac Low &; Ferraral (|l999h showed that metal- 
enriched material from stellar winds and supernova eject a 
is easily accelerated to velocities higher than the escape 
speed of the galaxy, especially in the case of low-mass sys- 
tems. This result, together with the fact that dwarf galaxies 
dominate the cosmic scenery in terms of number de nsity 
(|Springel et al.ll2005l : ISandage et ahlll985l : lMateolll998h . im- 
plies that dwarf galaxies play a crucial role in regulating the 
metallicity of the inter-galactic medium. 

The gravitational potential of NGC 2915 is dominated 
by its dark matter halo. Treating the density profile of the 
dark matter halo to be a pseudo iso-thermal sphere, the 
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Figure 6. Top panels: Position-velocity slices extracted from the model and Hi data cubes of NGC 2915. Each row represents a different 
position- velocity slice, with the slice locations and orientations in the H I total intensity map shown in the lower panel. 
Comments on columns: Each of columns 1-6 (from left to right) represent a different model cube of a thin, rotating gas disc. Columns 
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centre. Column 7 shows the Hi data cube position- velocity slices. 



escape velocity of the system is given by 



V e sc(R) 



2V C 2 1 + In 



Rn 



R 



(4) 



where V c is the asymptotic value of the circular rotation ve- 
locity at large radii, and R m ax is the maximum radius of 
the pseudo-isothermal sphere. Based on the rot ation curve 
derived for NGC 2915 by lElson et al.1 (|2010bh . we adopt 
V c - 86 km s" 1 and R max ~ 8-7 kpc. Using these values, an 
escape velocity radial profile was determined using Eqn. [4l 



For the sake of brevity, we simply mention that the escape 
velocities at all radii are > 175 km s _1 , much larger than 
the inferred expansion velocities of V exp ~ 30 km s _1 . Fur- 
thermore, the dark matter halo of NGC 2915 is undoubtedly 
larger than the radial extent of the gas disc, implying that 
these escape velocities are lower limits. The expanding gas 
in the centre of NGC 2915 therefore stands very little chance 
of escaping into the inter-galactic medium, and is consistent 
with a central outflow rather than a galactic blowout or 
wind. Assuming an outflow time-scale of ~ 100 Myr, the av- 
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erage mass redistribution rate is ~ 0.021 Mq yr 1 , roughly 
consistent with the estimated total star formation rate for 
the system, SFR to t ~ 0.028 M yr -1 (Section H2}. De- 
spite being unable to escape the gravitational potential of 
the galaxy, the redistributed gas will significantly affect the 
future star-forming properties of NGC 2915's stellar disc. 



escape velocities are a factor of ^ 6 larger than the estimated 
expansion velocities. Assuming a time-scale of ~ 100 Myr, 
the central gas is being expelled outwards at a rate that is 
roughly consistent with the estimated total star formation 
rate. The central H I kinematics of NGC 2915 should there- 
fore be thought of as forming a central outflow rather than 
a large-scale galactic blowout or wind. 



7 SUMMARY AND CONCLUSIONS 

We have used new deep, high-resolution Hi synthesis obser- 
vations of NGC 2915 to carry out a detailed study of its Hi 
kinematics which are known to be extremely complex within 
the immediate vicinity of the young stellar core. Sharp 
twists and kinks are seen in the iso-velocity contours of 
the Hi velocity field, and central Hi line profiles are known 
to be double-peaked and broad. We parameterised the Hi 
line profiles, checking which of them are better described 
as a double-component Gaussian as opposed to a single- 
component Gaussian. More than 20 per cent of the H I line 
profiles were found to be double-peaked, with many of these 
profiles co-located with the stellar disc of the galaxy. The 
average line splitting was measured to be AV = 59 km s _1 , 
which if interpreted in terms of an expanding gas component 
implies expansions velocities of ~ 30 km s _1 for roughly 
2.1 x 10 6 M of Hi. 

Ha emission is observed to spatially coincide with the 
identified double-peaked H I line profiles at inner radii. As- 
suming the deposition of mechanical energy from super- 
novae into the ISM, the energetics of a possibly radially- 
expanding inner H I component can be accounted for within 
~ 10 6 — 10 7 yr, the typical main sequence lifetimes of high- 
mass O or B type stars. The high-mass stellar disc near the 
centre of the galaxy is thus implicated as being the energis- 
ing source responsible for the anomalous H I kinematics. 

The main focus of this paper has been the construc- 
tion of full three-dimensional models of the NGC 2915 Hi 
data cube. These models were used to test various physi- 
cal scenarios which might be consistent with the system's 
observed Hi kinematics. Models based on purely circular 
dynamics are able to reproduce the observed kinematics of 
the outer Hi disc, yet fail to accurately reproduce those 
of the inner disc. Furthermore, the possibility of a high- 
velocity-dispersion inner portion of the H I disc, which may 
ordinarily be related to high levels of star formation activ- 
ity, is discredited. These variable-dispersion models are able 
to reproduce the large spreads in observed central velocities, 
yet cannot match the double-peaked structure of the H I line 
profiles. A final set of models constructed by adding radial 
expansion velocities 30 km s _1 at radii spanning the spa- 
tial extent of the stellar disc are able to reproduce the main 
kinematic features of the inner H I disc. The optimal model 
consists of a highly-inclined (i = 85°) inner Hi disc, without 
which the specific nature of the split H I profiles could not be 
matched. A highly-inclined inner disc could be the result of 
a dissipative merging process or perhaps a cold gas accretion 
event. Although direct evidence for such scenarios is lacking 
in the case of NGC 2915, they cannot be excluded as possi- 
ble causes of a highly-inclined inner disc. Due to the extreme 
dark matter halo properties of NGC 2915, our identified ex- 
panding central gas component stands essentially no chance 
of escaping into the inter-galactic medium. At all radii the 
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